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ANALYSIS OF INTEGRAL LIFT-FAN ENGINE DYNAMICS 
by John R. Szuch 
Lewis Research Center 

SUMMARY 

This report describes the hybrid computer simulation of an integral lift -fan engine. 
This engine is being considered for use in future VTOL airplanes. The integral lift -fan 
engine is, in principle, a conventional two -spool turbofan engine with a high bypass ratio. 
A contractor -proposed fuel control system and a simple model of the roll dynamics of a 
hovering VTOL airplane were added to the simulation. 

Both steady-state design and off -design data were generated using the simulation. 
Transient data were also generated to determine the effects of changes in the fuel control 
parameters, roll control parameters, and ambient temperature on the engine thrust and 
airplane roll angle response characteristics. 

The results of the hybrid computer studies indicated that the integral engine concept 
lends promise of achieving the response goals of a VTOL propulsion system. Although 
not optimized, the nominal fuel control gains and acceleration fuel schedule resulted in 
engine time constants less than 0. 2 second for standard -day thrust increments less than 
10 percent of the design thrust. The determining factor in achieving a 0. 20-second time 
constant for large thrust increments and/or hot -day conditions will be the tolerance of 
the turbines to transient overtemperatures. 

With the deceleration schedule adjusted to give "mirror -image" thrust responses, 
the nominal roll control gains resulted in unsatisfactory roll angle responses. A 50- 
percent increase in the roll rate feedback gain was required to achieve a satisfactory 
roll angle response. For a 10° roll angle demand, the increased feedback gain resulted 
in a 13 -percent overshoot in roll angle with 84 percent of the demand achieved in 1 sec- 
ond. For roll demands less than 5°, no overshoot was observed with 92 percent of the 
demand achieved in 1 second. 


INTRODUCTION 

Advances in turbofan engine technology and aircraft design have now made low -noise, 
high-performance VTOL and STOL intercity transports possible (refs. 1 and 2). NASA 



has been involved in both inhouse and contract studies of potential V/STOL propulsion sys- 
tems (refs. 3 to 5). A number of proposed VTOL airplanes employ turbine -driven lift 
fans. The lift fan and drive turbine may be an integral part of the gas generator or may 
be remotely located. Figure 1 illustrates an integral lift -fan propulsion system (ref. 6). 
The self-contained drive engine eliminates the ducting (hence, weight) associated with re- 
motely driven fan systems. The integral lift -fan engine is, in principle, a conventional 
two-spool turbofan engine with a high bypass ratio. The fan is shown with the core flow 
split off from the bypass flow. This geometry is attractive because it permits the hub 


Low-pressure 

compressor^ 

i 

i+ 



Main lift fan 


Hub fan 


r S ingle-stage high- 
\ pressure turbine 


(-Low-pressure 



'-Reverse-flow combustor 



uSS 




CS -56578 


Figure 1. - Integral lift-fan engine. 


section to operate at a lower pressure ratio than the main fan, thus permitting a smaller, 
lighter engine. The fan hub serves as a supercharger for the compressor. A reverse- 
flow combustor and a single-stage high-pressure turbine contribute to a short overall en- 
gine length. 

This report describes the analytical simulation of the integral lift -fan engine. An ex- 
isting general engine simulation program was used to implement the engine model on the 
hybrid (analog -digital) computer. A contractor -proposed fuel control system and a sim- 
ple model of a VTOL airplane (in the roll mode) were also added to the simulation. The 
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assumption of "mirror -image" thrust responses for accelerations and decelerations al- 
lowed rolling moments to be generated with only one engine being simulated. 

Both steady -state design and off -design data were generated using the simulation. 
Transient data were generated to determine (1) the effects of changes in the fuel control 
system on the engine thrust response characteristics, (2) the effects of changes in the 
ambient temperature on the engine response, (3) the effects of providing a turbine tem- 
perature limit in the fuel control on the engine response, (4) the acceleration and decel- 
eration schedules of fuel flow required for fast response and mirror -image thrust char- 
acteristics, (5) the effects of thrust level and thrust increment on the engine "time con- 
stant, " (6) the possible interactions between the airplane -engine control systems, (7) the 
effects of changes in the roll control system on the closed -loop response of the airplane, 
and (8) the effects of ambient temperature on the roll response characteristics. 

MATHEMATICAL MODEL 
Engine 

Figure 2 shows, in block diagram form, the model assumed for the integral lift -fan 
engine. As mentioned previously, the engine contains two spools or drive systems. A 
high-pressure turbine drives the compressor while a low-pressure turbine drives the 
supercharger - lift -fan combination. The supercharger (core) and fan (bypass) flows are 
physically split. A portion of the compressor discharge air is bled off for turbine cool- 
ing. Schematically, some of the cooling air reenters the cycle downstream of the high- 
pressure turbine (w^). The remainder of the cooling air (w^^) reenters the cycle 
downstream of the low pressure turbine. An overboard bleed flow (w ov ^) represents the 
flow lost from the cycle. 

Except for the bleed flows, all the compressor discharge air enters the combustor 
where it reacts with the injected fuel to generate the high energy gas needed to drive the 
turbines. The discharge from the low pressure turbine is ducted to a converging nozzle 
which accelerates the flow and gives axial thrust. The discharge from the fan is also 
ducted to a converging nozzle which provides about 80 percent of the total engine thrust. 

As indicated in figure 2, pressures and temperatures are computed in control vol- 
umes. In these volumes, storage of mass and energy occur (ref. 7). Pressure and tem- 
perature are assumed to be uniform in the volumes. The continuity equation, written for 
each volume, is 
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Figure 2. - Block diagram representation of integral lift-fan engine model. 

















(All symbols are defined in appendix A.) A summary of equations, as applied to specific 
components, is contained in appendix B. The energy equation, written for each volume, 
is 



When the results from equations (1) and (2) are used, the pressure in the volume can be 
computed from 


p _ RWT (3) 

V 

While the gas constant R is, in general, a function of the local fuel -air ratio, it was de- 
termined that the sensitivity of R to the fuel -air ratio could be neglected throughout the 
engine. 

In addition to the inter component volumes, flows must be computed in the combustor 
and core duct. For these cases, the pressure drop across the control volume is assumed 
to be due to both friction and momentum effects. The momentum equation (ref. 8) is 
written for each control volume as 

( p in- p om-**> + »i < 4 > 

For the generalized hybrid computer program, fans, compressors, and turbines are 
represented by overall performance maps. Available component data must be converted 
to the forms shown in figure 3. Inputs to the component maps are pressure ratio and a 
speed parameter. For fans and compressors, the map outputs are a corrected flow pa- 
rameter and adiabatic efficiency. For turbines, the outputs are a flow parameter and an 
enthalpy drop (work) parameter. A radial interpolation routine was used to generate the 
map outputs. For all components, the map outputs are used to compute flows, discharge 
temperatures, and torques. The computation of discharge temperature and torque, how- 
ever, requires knowledge of the "average” values of Cp and y in the component. 

These properties are assumed to be functions of temperature (the pressure dependence is 
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Figure 3. - Form of component performance data used in 
general engine simulation program. 



neglected). Therefore, they vary throughout the component. For this reason, an inter- 
polation constant /3 is adjusted for each component to match available cycle data; that is, 

T = ^T in + (1 - /3)T out (5) 




( 6 ) 

(7) 


r = 



( 8 ) 


where (f/a) is the local fuel -air ratio. 

After computing the fan, compressor, and turbine torques, the rotor speeds are de- 
termined by applying the principle of conservation of angular momentum; that is, 



(9) 


For the integral lift -fan engine, the supercharger characteristics did not require mapping. 
Because of the very low pressure ratio involved, the efficiency could be assumed constant 
and the pressure ratio approximated by a linear function of rotor speed. The volume be- 
tween the supercharger and compressor Vg j was neglected, and the compressor flow 
was used to compute the supercharger torque. For the integral lift -fan engine, flows 
must be computed for both the core nozzle and fan discharge nozzle. All nozzles are as- 
sumed to be the convergent type and the flow processes throughout the nozzles are as- 
sumed to be isentropic. For a specified inlet pressure, there exists for each nozzle a 
critical back pressure (ref. 7) given by 


P 


cr 



( 10 ) 
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If the back pressure on the nozzle is higher than the critical value, the flow is subsonic at 
the throat. For subsonic flow, the flow and thrust are given by 



If the back pressure is lower than the critical value, the flow at the throat is sonic or 
"choked. " For this case, the flow and thrust are given by 




y+l/y-1 


(13) 




+ A(P 


cr 


- P out^ 


(14) 


Because of the computing equipment limitations, it was decided to neglect the effects of 
variable y in the nozzle flow and thrust equations. A value of y was selected to best 
match available cycle data. The turbine cooling bleeds and overboard bleed flow were 
computed using equation (13). The low back pressures allowed the flows to be computed 
without testing the associated pressure ratio. 


Fuel Control 

Figure 4 shows a block diagram of a fuel control system proposed for the integral 
lift -fan engine. The function of the fuel control is to position a metering valve piston (x y ) 
to provide the proper amount of fuel to the combustor so as to maintain fan speed at a 
demanded value N f dem . Inputs to the fuel control system are fan speed N f , compres- 
sor speed N , compressor discharge pressure Pg, demanded fan speed Nj and 

maximum allowable compressor speed N c max . 

Transducer dynamics are included in the fuel control model. The two speed pickups 
are represented by 10-millisecond lags with gains of and Kg. The pressure trans- 
ducer is assumed to have a 30 -millisecond time constant. The metering valve is repre- 
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Figure 4 - Block diagram of contractor-proposed fuel control system for integral lift-fan engine. 


sented by a gain Ky. An actuator gain K AC relates the error in valve position to the 
rate of valve displacement k. An actuator gain Kg relates a speed error to the rate of 
valve displacement. Valve position is the integral of and is limited between 0. 2 and 
1. 0. Valve area is related to piston displacement by a square -law characteristic with the 
supplied fuel flow given by 

»F = < l5) 

Three MAX -MIN circuits are used to provide the logic needed to combine fast engine 
response with overspeed and overtemperature protection. The output of a MAX circuit is 
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equal to the largest of its inputs. Similarly, the output of a MIN circuit is equal to the 

smallest of its inputs. MIN circuit (A) is used to switch control to compressor speed if it 

approaches the maximum allowable speed N . MAX circuit (B) and MIN 

c y max 

circuit (C) determine whether valve position Xy is controlled by speed error or by the 
acceleration or deceleration fuel flow schedules. For example, in the initial stage of a 
deceleration, signals y ^ and y^ are negative. MAX circuit (B) selects y^ for its out- 
put which results in closed loop control of Xy with Xy dece j as the commanded position. 
Later in the transient y^ will exceed y^ and the valve position will be controlled by the 
speed error. 

The acceleration schedule, as shown, consists of a map of 


^max I - 



accel 


as a function of compressor speed N c with a fan speed -sensitive correction factor T2C. 

The deceleration schedule . is assumed to be equal to w /2. The measured 

pressure Pg m is used to multiply <? max and <p min to find the corresponding limits 
on fuel flow and valve position. 

To increase the system damping, hence stability, rate feedback is used to control the 
metering valve. While the block diagram indicates actual measurement of Xy with a 
feedback transfer function having a lag characteristic (i. e. , y 2 = K^y/K^s + 1)] x^, the 
actual feedback would probably be accomplished by passing a position signal Xy through 
a lead -lag network having the required frequency response characteristics. A summary 
of the equations describing the fuel control system is given in appendix B. 


Airplane Roll Control 

Figure 5 shows an illustration of the simple VTOL airplane configuration assumed 
for this study. The selected VTOL airplane weighed W^p at takeoff and used n E lift 
engines (n E /2 to a side). The airplane had a roll moment of inertia of I^p. 

The lift engines were all at an equal moment arm of r^ p and were used for both lift 
and roll control. For roll simulations, all engines were assumed to respond identically 
on each side with the opposite engines responding in a mirror -image fashion. This as- 
sumption allows roll -mode simulation with only one engine being simulated. While not al- 
ways achievable, the mirror -image characteristic would also be desirable from a fuel 
control implementation viewpoint. 

Figure 6 shows a block diagram of the airplane roll control system. The inputs to 
the roll control system are the demanded roll angle 6^ em , the measured roll angle e, 
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Figure 6. - Block diagram of integral engine VTOL airplane roll control. 


and the measured roll rate 6. The angular roll rate is fed back to add damping to the 
closed-loop system. The input to the fuel control dem is the sum of the throttle de- 
mand N f dem a and the roll angle error demand dgm To prevent over speeds 
and/or overtemperatures, the fuel control input was limited to dem < N f jqq where 
Nf jqq is the fan speed at the design thrust. For this study, roll maneuvers were con- 
ducted at the takeoff thrust and the fuel control input was limited to dgm > 0. 7 N f 
The limits would probably be temperature sensitive to account for changes in the ambient 
temperature. The fuel control generates the fuel flow w-p to the engine which, in turn, 
generates the thrust (lift) F. The n E /2 engine pairs generate a rolling angular accel- 
eration, which is given by 

.. n E r AP 

e = K e (f - F.) (16) 

I A Ti 
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where F. is the engine thrust which results in a steady nonrolling condition. The angu- 
lar roll rate and roll angle are determined by successive integrations of the angular ac- 
celeration. Both the roll angle and roll rate are fed back to the roll control system. 

The error between the demanded roll angle and the measured roll angle is multiplied by 
the gain Kj. The angular roll rate is multiplied by the gain Kg and subtracted from 
the amplified angle error. A temperature -sensitive gain is used to compensate for 
changes in the ambient temperature. Therefore, the roll error demand is given by 


N 


f,dem, e 1/ T 


2,des 


[ K l<»dem ' »> - K 2»] 


(17) 


The roll angle results in a lateral acceleration of the airplane given by 


n p F. 

X = K v -fL-i g c sin e 


W 


(18) 


AP 


For a hovering condition, the total thrust is equal to the airplane weight. 

The lateral velocity and displacement are determined by successive integrations of 
the lateral acceleration. A summary of the equations describing the airplane roll control 
model is given in appendix B. 


COMPUTER PROGRAM 

The integral lift -fan engine was simulated using an existing general engine program. 
The computer program allows simulation of either turbojet or turbofan engines while min- 
imizing the required manpower and programming effort. The use of the hybrid computer 
combines the accuracy associated with a digital computer with the continuous integration 
and user monitoring associated with the analog computer. 

The digital computer was used to (1) control the operation of the program, (2) per- 
form the logic required to set up the integral engine configuration, (3) perform the func- 
tion generation associated with the engine component models, and (4) perform most of the 
algebraic computations. The digital program was structured in subroutine form to min- 
imize core storage requirements and to allow convenient program modification if re- 
quired. The digital portion of the Lewis Research Center hybrid computer has 16 384 
words of core storage. Floating point (real) numbers are represented by two 16-bit 
words. Scaled fractions (less than 1. 0) are represented by a single 16-bit word. To 
minimize both the core requirements and computation time, the general engine program 
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was written using scaled -fraction digital computation. 

The analog portion of the hybrid computer provided the integration (with respect to 
time) associated with the system dynamics. The analog computer was used to compute 
stored weights, temperatures, pressures, duct flows, and rotor speeds. The analog 
computer was also used to implement the fuel control and airplane roll control models. 
Conventional strip-chart recorders and X-Y plotters were used to monitor the transient 
behavior of selected engine -airplane variables. 

For the integral lift -fan engine simulation, the analog portion of the program was 
time scaled (slowed down) by a factor of 100. This provided an acceptable ratio of digital 
cycle frequency to analog system frequencies and allowed the stepwise varying digital in- 
puts to the analog to be treated in a continuous manner. The combination of the 100:1 
time scaling and the digital cycle time of 43 milliseconds represents a 10:1 speedup in 
computation time when compared with all-digital simulation experiences. 


RESULTS 

The integral lift -fan engine (figs. 1 and 2) was simulated on the hybrid computer. 
Engine design point data, together with selected amplitude and time-scale factors, were 
used to generate input data required by the general engine simulation program. The en- 
gine design point corresponds to hot -day (305.4 K, 549. 7° R), maximum thrust condi- 
tions. Also put into the program were performance data for the eight engine component 
maps (two each for the fan, compressor, high-pressure, and low-pressure turbines). 


Steady-State Engine Operation 

The hybrid computer simulation of the integral engine was operated at the hot -day 
design point. The temperature interpolation constants (/3^, P c , *1 and ^g) were 
adjusted to achieve a steady-state cycle balance at the design point. 

The resultant steady-state operating point differed from the design point by less than 
0. 8 percent for system pressures, 0. 18 percent for system temperatures, 0. 3 percent 
for system flows, 0. 35 percent for rotor speeds, and 0. 65 percent for engine thrust. By 
varying the throttle demand a on the analog, the engine was throttled to the nominal 
takeoff condition (80 percent thrust) and selected engine variables were recorded. The 
ambient temperature Tg was then decreased to the standard-day temperature (288. 2 K, 
518. 7° R). The throttle demand was adjusted to give the same takeoff and maximum 
thrust values. Table I lists computed engine variables at the takeoff and maximum thrust 
conditions for both standard and hot days. The table also presents the end conditions for 
a number of transients run during this program . 
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TABLE I. - STEADY-STATE COMPUTED ENGINE DATA 


[All values specified as percent of design point. ] 



Hot day 

Standard day 

Takeoff 

Maximum thrust 

Takeoff 

Maximum thrust 

Thrust 

79.5 

99.4 

79.9 

99.6 

High-pressure 

92. 1 

99.9 

87.7 

95.6 

turbine inlet 





temperature 





Fan speed 

89.4 

100.3 

87.2 

97.6 

Compressor speed 

95.6 

100. 1 

93.0 

97.6 


Fuel Control - Engine Response 

To evaluate the fuel control system proposed for the integral lift -fan engine (fig. 4), 
accelerations from the takeoff condition to the design thrust condition were rim for both 
standard and hot days. It was decided that, for the range of rotor speeds involved, the 
acceleration schedule m could be approximated by a constant. Since the hybrid sim- 
ulation uses total pressures, the computed values of are approximately 4. 3 per- 

cent lower than the actual values (fuel control uses static pressures) at the design point. 

Effect of acceleration schedule . - Figure 7 shows the effect of changing the accelera - 
tion schedule <P mSiX on the responses of thrust and turbine inlet temperature to a step 
change in throttle position. Standard-day conditions were assumed. The figure also shows 
the tradeoff between the reduced time constant T(time to achieve 63 percent AF) and the 
increased peak turbine inlet temperature T^. The nominal value of (p max results in a 
time constant of 0. 285 second with a peak turbine inlet temperature equal to 99. 3 percent 
of the design value. The accepted response goal for a VTOL propulsion system is 
rsO.20 second (ref. 9). Figure 7(b) shows that turbine inlet temperature increases un- 
til the fuel control switches to closed -loop control of fan speed. The duration of turbine 
overtemperatures might prove to be as important a factor as peak temperatures. The 
tolerance of the integral engine to temperature excursions above the design value of T^ 
is undefined. The delay in achieving the final thrust value (see fig. 7(a)) is characteristic 
of the fuel control design and might be reduced by adding proportional action to the con- 
troller. This delay does not significantly degrade the overall performance of the engine. 

Figure 8 shows the effect of <P max on the thrust and turbine inlet temperature re- 
sponses, respectively, for the hot -day condition. Figure 8(a) indicates that a value of 
<p max = 04 ^ <p max nom is required to match the standard -day time constant of 

0. 285 second. The implementation of an ambient temperature -sensitive acceleration 
schedule might take the form shown in figure 4; that is, a temperature correction factor 
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High pressure turbine inlet Total thrust, percent of design 

temperature, percent of design 


Acceleration 

schedule, 

^max 


r \.QM Nominal 



(b) High-pressure turbine inlet temperature as function of time. 

Figure 7. - Effect of acceleration schedule <p m ax on engine response. Accelerations from takeoff to 
design thrust; standard-day, sea level static conditions; fuel control gains nominal. 



Acceleration schedule, 
'"max 



(b) High-pressure turbine inlet temperature as function of time. 

Figure 8. - Effect of acceleration schedule <p m3X on engine response. Accelerations from takeoff to de- 
design thrust; hot-day, sea level static conditions; fuel control gains nominal. 


T2C would be computed from the measured fan speed and used to modify the acceleration 
schedule <p (N ). Figure 8(b) shows a peak turbine inlet temperature equal to 1. 051 
times the design value resulting from the increased ambient temperature and acceleration 
limit. 

To evaluate the sensitivity of the engine response to changes in the fuel control gains, 
standard-day accelerations from the takeoff condition to the design thrust condition were 
run. The nominal value of <p max was used for these runs. 

Effect of valve actuator gain . - The valve actuator gain Ky (see fig. 4) was varied 
from 0. 2 to 4. 0 times the nominal value. This gain affects both the initial phase of the 
acceleration (closed loop on valve position) and the latter phase (closed loop on fan speed). 
Figures 9(a) and (b) show the resulting standard -day responses of thrust and turbine inlet 
temperature, respectively. The lowest gain results in a sluggish response (t = 0. 52 sec- 
ond) with little overshoot in temperature. The highest gain resulted in a reduction in the 
time constant to r = 0. 235 second at the cost of slightly increased peak temperature 
(1. 007 times design value) and a more oscillatory response. 

Effect of valve rate feedback gain. - To increase the damping of the fuel control sys- 
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Metering valve 



(b) High-pressure turbine inlet temperature as function of time. 

Figure 9. - Effect of fuel valve actuator gain K v on engine response. Accelerations from takeoff to 
design thrust; standard-day, sea level static conditions; acceleration schedule, p max ■ nominal; all 
other fuel control gains nominal. 


tem, the rate feedback transfer function K^y/(r^s + 1) is used to feed back the valve 
rate to the controller. To evaluate its effect on the engine response, the gain K^y was 
varied from 0. 1 to 2. 0 times the nominal value for the standard -day accelerations. Fig- 
ures 10(a) and (b) show the thrust and temperature responses for this case. They also 
show that the rate feedback only affects the closed -loop response after the fuel control 
switches from the acceleration schedule to the closed -loop control of fan speed. Decreas- 
ing K^y results in slightly increased peak turbine inlet temperature while delaying the 
switch to fan speed control. As expected, decreased K^y also results in a more oscil- 
latory response. The nominal value results in a satisfactory response although K^y 
might be decreased to hasten the achievement of maximum thrust at the cost of decreased 
damping. 

Effect of temperature limiting . - Figures 8(a) and (b) indicate that significant turbine 
overtemperatures result when fast, hot -day accelerations to maximum thrust are run. 
Most turbojet and turbofan engines incorporate overtemperature limits in their control 
systems. To evaluate the response degradation associated with temperature limiting, the 
circuit shown in figure 11 was added to the fuel control system (see fig. 4). MIN circiut 
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High-pressure turbine inlet Total thrust, 

temperature, percent of design percent of design 


Rate feedback gain 
for fuel valve, 




(b) High-pressure turbine inlet temperature as function of time. 

Figure 10. - Effect of fuel valve rate feedback gain K^y on engine response. Accelerations from take- 
off to design thrust; standard-day, sea level static conditions,- acceleration schedule, (o max • nominal; 
all other fuel control gains nominal. 



i 


Figure 1L - Block diagram of temperature override in fuel control. 








High-pressure turbine inlet Total thrust, High-pressure turbine inlet Total thrust, 

temperature, percent of design percent of design temperature, percent of design percent of design 


Temperature ratio, 
t 4, max 
^4,des 



(c) Effect of gain K 6 on thrust response for temperature ratio max /T 4 des of 1. 03. 



Time, sec 


Id) Effect of gain on high-pressure turbine inlet temperature response for temperature ratio 
^4,max^4,des 

Figure 12. - Effect of temperature T 4 limit on engine response. Accelerations from takeoff to design 
thrust; hot-day, sea level static conditions; acceleration schedule, (o max - 1.044 nominal; fuel 
control gains nominal. 



(D) was used with a gain Kg controlling the switch to closed -loop control of the inlet tem- 
perature T 4 if the maximum allowable temperature T. is approached. Fig- 

ures 12(a) and (b) show the effect of the temperature limit on the hot -day responses of 
thrust and turbine inlet temperature, respectively. The acceleration schedule <p max 

was equal to 1. 044 times the standard -day value for these runs. The gain K* was set to 
-3 b 

5. 055x10 for the accelerations shown in figures 12(a) and (b). With the maximum allow- 
able temperature T 4 max set to 1. 10 times the design point value, the engine responses 
correspond to those obtained with no temperature limiting (see fig. 8). As T d ap- 

proaches the design point value, however, the thrust response becomes more sluggish. 

A value of T 4 max equal to 1. 03 T 4 deg resulted in a 0. 4-second time constant with 
the greater portion of the acceleration controlled by T 4 . With T, equal to 

T 4 des , the time constant increased to 0. 735 second. Figures 12(c) and (d) show the ef- 
fect of the gain Kg on the thrust and turbine temperature, respectively. Temperature 
T 4 max was set to 1.03 T 4 deg for these runs. Since increasing K g delays the switch 
to temperature control (in addition to increasing the control gain), some reduction in the 
time constant is observed. With Kg equal to 1. OxlO -2 the time constant was 0. 34 sec- 
ond and the peak temperature was very nearly equal to the specified maximum. 

While the fuel control gains were not optimized for this portion of the study, it is ap- 
parent that significant increases in the engine time constant can result from temperature 
limiting. The need for such limiting will be based on the tolerance of the turbines to the 
peak temperatures associated with rapid accelerations. 

Effect of deceleration schedule . - Having determined the acceleration characteristics 
of the integral engine and fuel control, a study was undertaken to determine the mirror - 
imageness of the engine’s thrust responses. 

While the engine manufacturer had assumed a deceleration schedule w . equal to 

nmn 

one-half the acceleration schedule <P max , the significance of this assumption was un- 
known. For comparison with available acceleration data, decelerations were run from 
the takeoff condition (0. 8 F deg ) to the 60- per cent thrust condition. Standard-day condi- 
tions were assumed with <p max equal to the nominal value and ^ min /<p max varied from 
0. 5 to 0. 7. Figure 13 shows the resultant thrust responses. A ratio of 0. 65 results in 
the desired mirror -image characteristic. Similar decelerations were run for the hot -day 
condition with (p max equal to 1.044 times the nominal value. Figure 14 shows the re- 
sultant thrust responses. The ratio of 0. 65 proved satisfactory for the hot -day condition, 
also. Figure 15 illustrates the effect of <P m i n /<P max on the deceleration time constant 
for both standard -day and hot -day conditions. 

Effect of thrust level and thrust increment . - The response time of the integral en- 
gine is a function not only of the fuel control parameters but also of the thrust level and 
thrust increment. Standard -day accelerations from initial thrusts, ranging from the 45- 
percent level to the takeoff thrust (80 percent), were run using the nominal fuel control. 
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Figure 14. - Effect of deceleration schedule <p m - m on engine thrust response. Decelerations from 
takeoff to 60-percent thrust; hot-day, sea level static conditions; fuel control gains nominal. 



Figure 15. - Effect of deceleration schedule g> m j n on engine time 
constant. Decelerations from takeoff to 60-percent thrust; fuel 
control gains nominal; sea level static conditions. 
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Figure 16. - Effect of thrust increment on engine response. Accelerations from landing thrust; standard-day, 
sea level static conditions; acceleration schedule, y> max = nominal; fuel control gains nominal. 


The assumption of <P m i n /V max = 0 . 50 was used to allow comparisons with contractor - 
generated data. Decelerations were also run from thrust levels ranging from the design 
thrust to 65-percent thrust. Figures 16(a) and (b) show the responses of thrust and tur- 
bine inlet temperature, respectively, for the accelerations from the landing condition (52- 
percent thrust). The time constants for all accelerations and decelerations are plotted in 
figure 17 as functions of the initial thrust level and the thrust increment AF. Time con- 
stants as low as 0.095 second (for the 10 -percent F^gg deceleration from the design 
thrust) and as high as 0. 53 second (for the 50-percent F^ acceleration from the 55- 
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Normalized 



Figure 17. - Effect of initial thrust and thrust increment on engine 
time constant. Standard-day, sea level static conditions; accelera- 
tion schedule, <p mx * nominal; deceleration schedule, p m jn = 

0. 5?>max : fuel contr < )l 9 ains nominal. 


percent thrust level) were observed using the simulation. Figure 17 indicates that the de- 
sired response goal (r < 0. 20 sec) is achieved for thrust increments less than 10 percent 
of the design thrust. Most accelerations and decelerations during flight would involve 
thrust increments of this magnitude. 


Roll Control - Airplane Response 

10° Roll maneuver . - To aid in ascertaining the response interactions of the integral 
engine fuel control with a VTOL airplane, the aircraft roll control (fig. 6) was added to 
the hybrid simulation. A standard-day hovering condition at takeoff thrust was assumed. 
To allow possible comparisons with contractor data, the ratio (P m ^ n / <p max was se * *° 

0. 5. With the engine operating at the takeoff condition (80-percent thrust), a 10° 

(0. 1745 rad) roll angle demand was input to the roll control system. Figures 18(a) and 
(b) show the resulting responses of roll angle Q and engine thrust, respectively. All 
controller gains were nominal with <P max equal to the nominal value. A peak roll angle 
of 11.9° (0. 2077 rad) was observed with a roll angle of 8.9° (0. 1553 rad) achieved in 
1 second. A goal of reaching the demanded roll angle in 1 second with less than 10- 
percent overshoot is desirable for a VTOL airplane. Figure 18(a) indicates that this goal 
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might be achieved with some modification in the airplane and/or control system. Some 
modification to the roll control system might also be required to eliminate the 0.25- 
second deadtime observed in the roll angle response. Figure 19 shows strip-chart re- 
cordings of selected engine and roll control variables for the standard -day roll maneuver. 

Effect of deceleration schedule . - The same roll maneuver was run with the value of 
^min^^max = This value results in the assumed mirror -image thrust responses 

for accelerations and decelerations from the takeoff condition. Figures 18(a) and (b) 
show the effect of the increased <p min on the roll angle and thrust responses, respec- 
tively. The increased fuel flow results in a slower deceleration in response to the roll 
angle overshoot. A peak roll angle of 14. 6° (0. 2548 rad) results with a significant de- 
crease in the system damping. Figures 18(a) and (b) indicate that the selection of control 
gains must account for any modifications made in the fuel control system. 


i.o 








Figure 19. - Response of selected engine-airplane variables to 10° roll angle demand. Hover- 
ing at takeoff thrust; standard-day, sea level static conditions; acceleration schedule, 

(£> r()ax = nominal; deceleration schedule, ^ m j n = 0. 5<P max ; fuel control and roll control 
gains nominal. 
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Figure 20. - Effect of roll rate feedback gain K 2 on airplane response to 10° roll angle demand. Hover- 
ing at takeoff thrust; standard-day, sea level static conditions; acceleration schedule, y> ma)! ■ nominal; 
deceleration schedule, (o m j n * 0. 65(0 max ; all other fuel control and roll control gains nominal. 


Effect of roll rate feedback gain . - To increase the closed-loop system damping with 
the desired ^ min /<P max = 0.65, the roll rate feedback gain K 2 was varied from the 
nominal value to twice the nominal value. Figure 20 shows the resulting roll angle re- 
sponses for the standard -day , 10° (0. 1745 rad) roll maneuver. A 50-percent increase in 
the nominal roll rate feedback gain results in a peak roll angle of 11.3° (0. 1972 rad) with 
a roll angle of 8.4° (0. 1466 rad), achieved in 1 second. 

Effect of roll angle demand . - As in the case of throttle demands, the response char- 
acteristics of the roll control system are affected by the size of the input demand. Fig- 
ure 21 shows the response of roll angle to roll angle demands ranging from 3. 75° to 10° 
(0. 0654 to 0. 1745 rad) for standard-day, takeoff conditions. The increased roll rate 
feedback gain was used for these runs. The effective damping of the closed -loop system 
decreases as the roll angle demand increases. There is virtually no overshoot in the 
roll angle for input demands less than 5° (0.0873 rad). The roll angle, achieved in 1 sec- 
ond, varied from 84 percent (10° demand) to 92 percent (less than 5° demand) of the de- 
manded roll angle. These results are encouraging since a typical hovering maneuver 
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Time, sec 


Figure 21 . - Effect of roll angle demand on airplane response. Hovering at takeoff thrust; standard-day, 
sea level static conditions; acceleration schedule, (P max = nominal; deceleration schedule, g> m j n • 

0.65^max : 9 ain ' K 2 ' 15 nominal ; al1 other fuel control and roll control gains nominal. 

might involve roll angle demands of ±3° (±0. 0524 rad). 

Effect of ambient temperature . - As mentioned previously, a determining factor in 
the achievement of fast response will be the tolerance of the integral engine to the poten- 
tial turbine overtemperatures associated with hot -day accelerations of the engine. Since 
no contractor data for hot -day operation were available, the 10° (0. 1745 rad) roll maneu- 
ver was simulated at the hot -day takeoff condition with the increased roll rate feedback 
gain. The increased ambient temperature was reflected in (1) increased throttle demand 
Nf dem to achieve the takeoff thrust (80 -percent thrust), (2) increased roll error de- 
mand dem , (3) increased limits on the fuel control input demand dem , and (4) 
increased acceleration schedule <p max - Figures 22(a) and (b) show a comparison of the 
standard -day and hot -day responses of roll angle and turbine inlet temperature, respec- 
tively. Figure 22(a) shows that virtually identical responses of 6 are achieved by the 
temperature compensations previously described. Figure 22(b) shows the expected in- 
crease in turbine temperature due to the hot -day condition. In addition to the peak tem- 
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(b) High-pressure turbine inlet temperature as function of time. 

Figure 22. - Effect of ambient temperature on engine-airplane response to 10° roll angle demand. Hover- 
ing at takeoff thrust; sea level static conditions; deceleration schedule, <p m - m - 0. 65p max ; gain, 

«2 ■ 1. 5 nominal; all other fuel control and roll control gains nominal. 



perature (1. 054 times the design point temperature), the rapid changes in turbine tem- 
perature associated with the 10° roll maneuver might have an undesirable effect on the 
engine. 


SUMMARY OF RESULTS 

An integral lift -fan engine, proposed for use in a VTOL airplane, was simulated on 
the hybrid computer using a general engine simulation program. A fuel control system 
and a VTOL airplane roll control system were added to the engine simulation. Sea level 
static conditions were assumed. The results of the hybrid computer studies indicated the 
following: 

1. The nominal acceleration schedule results in a 0.285-second time constant for a 
standard-day takeoff to design thrust acceleration. The desired response goal of 

t < 0. 20 second was achieved for thrust increments less than 10 percent of the design 
thrust. 

2. Increasing the standard-day acceleration schedule 4.4 percent results in a 

0. 285- second time constant for the hot -day takeoff to design thrust acceleration; the in- 
creased ambient temperature and fuel flow results in a peak turbine inlet temperature 
5. 4 percent higher than the design point value. 

3. Although not optimized, the nominal fuel control gains result in satisfactory thrust 
response characteristics. An increased fuel valve actuator gain resulted in slight reduc- 
tions in the time constant at the cost of increased peak turbine temperature. 

4. The use of a temperature limit in the fuel control can seriously degrade the hot- 
day acceleration to design thrust. Limiting the turbine inlet temperature to the design 
point value increased the time constant from 0. 285 to 0. 735 second. 

5. The use of the contractor's assumed relationship between the deceleration and ac- 
celeration schedules (<P min = 0. 5 <P max ) does not result in the desired mirror -image 
thrust responses. A ratio of ^ m i n /V max = 0- 65 proved to be satisfactory for both 
standard -day and hot -day conditions. 

6. The use of ^ m j n /V max = 0. 65 results in an unsatisfactory roll angle response 
with the nominal roll control gains. For a standard- day 10° roll demand, a 46- percent 
overshoot in roll angle was observed indicating significant interaction between the fuel 
control and roll control systems. A 50- percent increase in the nominal roll rate feed- 
back gain with <P min = 0.65 <P max resulted in a satisfactory roll angle response. For 
a 10° roll angle demand, a 13-percent overshoot was observed with 84 percent of the 
demand achieved in 1 second. For roll demands less than 5°, no overshoot was ob- 
served with 92 percent of the demand achieved in 1 second. 



7. The use of temperature-sensitive throttle gain, roll control forward gain, accel- 
eration schedule, and fuel control input limits results in virtually identical roll -angle re- 
sponses for both standard -day and hot -day conditions. 

The results of this study indicate that the integral engine concept lends promise of 
achieving the response goals of a VTOL propulsion system. The determining factor in 
the achievement of the desired 0. 20 -second time constant for large thrust increments will 
be the tolerance of the turbines to overtemperatures during accelerations to design thrust. 
Once this tolerance is known, acceleration schedules, controller gains, etc. can be opti- 
mized to accomplish the necessary reduction in response time. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 19, 1972, 

764-72. 
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APPENDIX A 


SYMBOLS 

2 2 

A effective cross-sectional area, cm (in. ) 
c specific heat at constant pressure, J/kg-K (Btu/lbm-°R) 

c y specific heat at constant volume, J/kg-K (Btu/lbm-°R) 

cnp compressor speed parameter 

ctr compressor ideal temperature rise parameter 

d slope of supercharger performance map, N/cm (lbf/in. ) 
dt differential time, sec 

e intercept of supercharger performance map, N/cm (lbf/in. ) 

F thrust, N (lbf) 

^ functional relationship 

AF thrust increment, N (lbf) 
f/a fuel -air ratio 

fcnp general fan or compressor map speed parameter 
fnp fan speed parameter 

fpc compressor map flow parameter, kg/sec (lbm/sec) 
fpf fan map flow parameter, kg/sec (lbm/sec) 

fpfe general fan or compressor map flow parameter, kg/sec (lbm/sec) 

fpt high-pressure turbine map flow parameter, kg-K-cm /N-rpm-sec (lbm- R-in. / 
lbf-rpm-sec) 

2 o 2 

fptl high-pressure turbine map flow parameter, kg-K-cm /N-fpm -sec (lbm- R-in. / 
lbf-rpm-sec) 

2 o 2 

fpt 2 low-pressure turbine map flow parameter, kg-K-cm /N-rpm-sec (lbm- R-in. / 
lbf-rpm-sec) 

ftr fan ideal temperature rise parameter 

2 2 

g c gravitational conversion factor, 100 cm-kg/N-sec (386. 3 lbm-in. /lbf-sec ) 

HVF heating value of fuel, JAg (Btu/lbm) 

Ah turbine enthalpy drop, J/kg (Btu/lbm) 
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general turbine map enthalpy drop parameter, J /kg -K ^ 2 -r pm (Btu/lbm -°R^ 2 - 
rpm) 

high-pressure turbine map enthalpy drop parameter, J/kg-K*/ 2 -rpm (Btu/lbm- 
°R*/ 2 -rpm) 

low-pressure turbine map enthalpy drop parameter, J/kg-K^ 2 -rpm (Btu/lbm- 
°R 1 / 2 -rpm) 

2 2 

polar moment of inertia, N-cm -sec (in. -lbf-sec ) 

work conversion factor, 100 N-cm/J (9339. 6 in. -lbf/Btu) 

valve position controller gain 

rate feedback gain for fuel valve 

metering valve flow coefficient 

metering valve actuator gain 

distance conversion factor, 0.0833 ft/in. (1.0) 

roll controller forward gain, rpm/deg / 

roll controller rate feedback gain, rpm -sec/deg 

speed controller forward gain 

fan speed controller gain 

compressor overspeed governor gain 

turbine overtemperature governor gain 

angle conversion factor, 57.32° 

torque, N-cm (in. -lbf) 

differential torque, N-cm (lbf -in. ) 

length, cm (in. ) 

number of flows fed by control volume 
rotational speed, rpm 
number of flows feeding control volume 
number of engines 

2 

total pressure, N/cm (psia) 
static pressure, N/cm (psia) 
pressure ratio 
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prc compressor map pressure ratio 
prf fan map pressure ratio 

prtl high-pressure turbine map pressure ratio 
prt2 low-pressure turbine map pressure ratio 
R gas constant, N-cm/kg-K (in. -lbf/lbm-°R) 

0t duct pressure loss coefficient, N-sec Ag -cm (lbf-sec /lbm -in. ) 
r AP r °U momen t arm, in. (cm) 

s Laplace operator, sec"^ 

T total temperature, K (°R) 

T2C temperature -sensitive corrector factor on acceleration schedule 

t time, sec 

tnp general turbine map speed parameter, rpm/K 1//2 (rpm/^R 1 ^ 2 ) 

tlnp high-pressure turbine map speed parameter, rpm/K^ 2 (rpm/°R*/ 2 ) 

t2np low-pressure turbine map speed parameter, rpm/K 1 / 2 (rpm/°R 1//2 ) 

tr general fan or compressor ideal temperature rise parameter 

trsc supercharger ideal temperature rise parameter 
3 3 

V volume, cm (in. ) 

W mass, kg (lbm) 

w mass flow, kg/sec (lbm/sec) 

X aircraft lateral displacement, cm (ft) 

Xy valve displacement, cm (in. ) 

y i general fuel control variable (appropriate units) 

a relative throttle position 

/3 temperature interpolation constant 

y specific heat ratio 

r] efficiency 

e roll angle, deg (rad) 

r time constant, sec 

Tj fuel valve rate feedback time constant, sec 

2 

<p fuel flow parameter, kg-cm /N-sec (lbm/sec/psi) 



Subscripts: 

A air 

AP airplane 

a ambient 

accel acceleration 

b burner 

bll high-pressure turbine cooling bleed 

bl2 low-pressure turbine cooling bleed 

c compressor 

cd core duct 

cr critical 

decel deceleration 

dem demand 

des design 

F fuel 

f fan 

1 initial condition 

m measured 

max maximum 

min minimum 

nom nominal 

ovb overboard bleed 

SC supercharger 

std standard day 

tl high-pressure turbine 

t2 low-pressure turbine 

tot total engine 

a throttle demand 

9 roll demand 

2 fan inlet 
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2. 1 supercharger -compressor intercomponent volume 

2. 5 fan duct volume 

2. 8 fan nozzle volume 

3 compressor -combustor intercomponent volume 

4 combustor volume 

5 high-pressure turbine/low -pres sure turbine intercomponent volume 

5. 1 one -half of core duct volume (upstream) 

6 one -half of core duct volume (downstream) 

8 core nozzle throat 

70 70-percent maximum thrust 

100 100 -percent maximum thrust 

Superscripts: 

_ average 

f 

inlet of control volume 
time derivative 



APPENDIX B 


SUMMARY OF EQUATIONS 

Engine: 


prf = 2 ‘ 5 

P 2 

(Bl) 

V N f des 

fnp - 1 1>aes 

V T 2 //T 2,des 

(B2) 

fpf = (prf, fnp) 

(B3) 

A f P f ( P 2/ P std) 

w 2. 5 - , 

V^Std 

(B4) 

77f = ^2 f(P rf > f np) 

(B5) 

Vf = , 0) 

(B6) 

c p f =^6 (f f 0) 

(B7) 

= ^f T 2 + ^ ' %' T 2. 5 

(B8) 

/P 

ftr -( 2 ' 5 ) - 1.0 

\ p 2 / 

(B9) 

T 2.5=[(-V 10 1 T 2 

A^f/ J 

(BIO) 

“ C p f (T 2.5 " T 2 )w 2.5 

*' 

(Bll) 
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Cp 2 5 = ' r 6 (T 2.5’ 0) 

(B12) 

P 2 i = d(fnp) + e 

(B13) 

/p \0. 286 

trsc 2 - -1.0 

(B 14) 

t 2 i-[(— U1-o1t 2 

LV^sc J 

(B15) 

T T C Psc' T2 ' > ' T ^ W 3 

80 " N f 

(B16) 

prc = P 3 /P 2 j 

(B17) 

cnp - N c /N c,d es 

V T 2.1 /T 2.1,des 

(B18) 

fpc = ^i )C (P rc ’Cnp) 

(B 19) 

_ fpc ( p 2.l/ p std 

V T 2. l^std 

(B20) 

>J C = ^ c (prc, cnp) 

(B21) 

v c = 

(B22) 

c p c = ^c-°> 

(B23) 

? c = e c T 2. X + <1 - f! c )T 3 

(B24) 
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T 3 = 



+ 1.0 


2.1 


L c * 


™ c p c (T 3- T 2.1^ 


N 


% - *6< T 3’°> 


tlnp = 


N. 


ft 

prtl = P 5 /P 4 
fptl - 3 ^ ^(prtl,tlnp) 

. fetlPjN 


W 5 ~ 


hptl = ^ 4 )tl (prtl,tlnp) 
hptlNc 


Ah tl = 


c n = 

P U 6 


1000 


T trf~ 


4J 


T tl * % T 4 + " % T 5> 


30J 


L u = 


^5 


N„ 


(B2t 

(B27) 

(B28) 

(B29) 

(B30) 

(B31) 

(B32) 

(B33) 

(B34) 

(B35) 

(B36) 

(B37) 
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Ah+i 

T’-T, 1.0- 

5 4 l (wi 

(B38) 


(B39) 

(L) . !? 
W 4 w 4 

(B40) 

t2np = Nf 

(B41) 

prt2 = Pg i/ P g 

(B42) 

fpt2 = ^ t 2(prt2,t2np) 

(B43) 

fpt2P cN, 

W 5. 1 T 

*5 

(B44) 

hpt2 = ^ t 2(P rt 2,t2np) 

(B45) 

Ak h P t2N fV f I 

^^■0 

^ 1000 

(B46) 

c P.2^fMD 5 ’ 

(B47) 

*t2 = /3 t2 T 5 + ^ " /3 t2 )T 5. 1 

(B48) 

V Ah t2 W 5.1 
L t2 = JL -^ 

(B49) 


41 



(B50) 



y 2.5 " ^ T 2.5’ 0) 


° p 2 5 ~ '^ (T2 - 5 ’ 0) 


S.s'VsA 2 - 5 

y 3 = ^5^3’°* 


c p 3 = ^ T 3’ 0) 



y 4 = ^5 



(B51) 


(B52) 


(B53) 

(B54) 

(B55) 

(B56) 

(B57) 

(B58) 

(B59) 

(B60) 
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(B61) 

= % 

(B62) 

> 


T b - % t 3 + d - V T 4 

(B63) 

% = *6<V°> 

(B64) 

'• ' ''•[’••©,] 

(B65) 

C P g= 

(B66) 

c 

c -2® 

(B67) 

5 y 5 


r 5-t = # l[ T 5.1’(9 5 t 

(B68) 

%.r f2 [ T5 -'’(a) 51 _ 

(B69) 

c 

c = **■' 

V5 - 1 r 5.1 

(B70) 


(B71) 

V '•['•'©„] 

(B72) 
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Cy 6 ' r, 


P_6 

6 


(B73) 


T 3 = 


N,= ^X t(Lt2 ‘ Lf " LscWt+Nf .‘ 




-f 


W 3 = L (w 3 - w bll - w bl2 - w ovb - w 4>® + W 3,i 


(B74) 


(B75) 


(B76) 


f 

1 W A % 


*\ 


T 3 [»3 + fr 3 - Ufrt * *M1 + »bl2 + *ort>>] 


ldt + T 3i (B77) 


p R A W 3 T 3 

y 

V 3 


W, 


=^ t( -i 4 + w F -i 5 )dt+ w 4]i 


(B78) 


(B79) 




f( c p b T »'"4 + ”b HVFi r) 


- T 4 [w 4 + W F + (y 4 - 1)W 5 U + T 4}i (B80) 


P 4 = 


R A W 4 T 4 


(B81) 


W c 


‘ jC (w 5 + w bll - w 5. l Wt + W 5,i 


(B82) 
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T 5 = 


/ 

•'O 


1 J( C P 5 T 5-5 - c p 3 T 3-bll) 


W c 


- T 5 [*5 + *bll + ^5 - 1)W 5. l]f dt + T 5,i ( B83 ^ 


RaW.t,. 
p - A 5 5 


(B84) 


W c 


l’f, <*S.l + *btt-*6 )dt+W 5.1,i 


-I- 


r ( Cp 5 l^- 1 W 5. 1 + C p 3 T 3 W bl2 


•) 


'5. 1 


(B85) 


' T 5. l[*5. 1 + w bl2 + fr 5 . 1 - 1)W 6] 

RaW- ,t 


ydt + T 


5. 1, i 


5. 1 


A 5. 15. 1 

Vl 


W 6 = /' t(,i '6- i 8 Wt+W 6,i 


(B86) 


(B87) 


(B88) 



T 6 [^6 + ^6.- 1)W 8]f dt + T 6,i 


(B89) 


P 6 = 


R A W 6 T 6 


(B90) 


W, 


!.5=/ ( *2. 


5 - W 2.8 )dl+W 2.5,i 


(B91) 
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(B92) 



w F = output from fuel control 



(B93) 

(B94) 

(B95) 

(B96) 

(B97) 

(B98) 

(B99) 



(B100) 


w 


2.8 





(B101) 
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^ 7 (r) = 


y+l/y-1 


(B102) 



y + l. 




^(y,pr) = &j(y) if pr < 0.53685 

= (pr) 1/y ^ - y — jl - p r y ~ 1/y ] if pr >0.53685 


F S = '8 V». T V'* f « -H + V P ., S - P .' 




F 2.8 ~ w 2. 8 V c p 2 5 T 2.5 ~ ^9 f y 2.5 + A 2. 8 (P s, 2. 8 " P a ) 

c \ 2.5 


F tot ~ F 2. 8 + F 8 


P g 8 = 0. 53685 Pg if P a /P 6 =£ 0. 53685^ 

= P a if P a /Pg > 0.53685 

P s, 2.8 = 0 - 53685 P 2.5 if V P 2.5 s0 ‘ 53685 | 
= P a if P a /P 2- 5 > 0. 53685 

•?g(y,pr) = Vl - pr y-1//y 


Fuel control: 




= MIN< 


l K 4( N f,dem- N f, m ^ 
l K 5^c,max- N c,m) 


(B103) 

(B104) 

(B105) 

(B106) 

(B107) 

(B108) 

(B109) 

(B110) 
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for y 2 ^ 0. 8 


(BUI) 




/' 


(K 


AV x v 


- y 2 )dt 


y 3 = Kgfrj - y 2 ) 




*F - K mv4 
^max _ ^10^c,rn’^f, 
^min " ^ll^max^ 


(B112) 

(B113) 

(B114) 

(B115) 

(B116) 

(B117) 

(B118) 

(B119) 

(B120) 

(B121) 
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Airplane roll control: 


N M e m =\de m ,^ N f,dem, 9 ,OT \k> 0 S \de m * N l ,100 = N f,des 


(B122) 


*° r05 “ S1 -° 

(B123) 

N f,dem,e-^ T [ K l (g dem e) 2] 

’ " i 2,des 

(B124) 

e = K " ErAP (F - F t ) 
6 *AP 

(B125) 

9 = P 6 dt 
J 0 

(B126) 

e= r e dt 
-'O 

(B127) 

n E F.g sin 0 
v r ^ 1 c 

(B128) 

A ur 

W AP 


rt .. 

X = / X dt 

Jo 

(B129) 

X = f* X dt 
./n 

(B130) 
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